HMC algorithm for two-flavour lattice QCD: Schwarz-preconditioning with a one-dimensional domain decomposition
نویسنده
چکیده
We study a variant of the Schwarz-preconditioned HMC algorithm. In contrast to the original proposal of Lüscher, we apply the domain decomposition in one lattice direction only. This is sufficient to reduce the condition number of the fermion matrix restricted to the domains compared with the full fermion matrix. For the same linear extension of the domain, less links reside on the boundaries of the domains. Therefore it becomes e.g. practical to iterate the decomposition. We perform numerical tests for two degenerate flavours of Wilson fermions. The standard Wilson gauge action at β = 5.6 is used. The performance of our implementation is compared with other recent studies using various types of preconditioning.
منابع مشابه
Schwarz-preconditioned HMC algorithm for two-flavour lattice QCD
The combination of a non-overlapping Schwarz preconditioner and the Hybrid Monte Carlo (HMC) algorithm is shown to yield an efficient simulation algorithm for two-flavour lattice QCD with Wilson quarks. Extensive tests are performed, on lattices of size up to 32×24, with lattice spacings a ≃ 0.08 fm and at bare current-quark masses as low as 21MeV.
متن کاملSchwarz-preconditioned HMC algorithm for two-flavor lattice QCD
The combination of a non-overlapping Schwarz preconditioner and the Hybrid Monte Carlo (HMC) algorithm is shown to yield an efficient simulation algorithm for two-flavour lattice QCD with Wilson quarks. Extensive tests are performed, on lattices of size up to 32×24, with lattice spacings a ≃ 0.08 fm and at bare current-quark masses as low as 21MeV.
متن کاملHMC algorithm with multiple time scale integration and mass preconditioning
We present a variant of the HMC algorithm with mass preconditioning (Hasenbusch acceleration) and multiple time scale integration. We have tested this variant for standard Wilson fermions at β = 5.6 and at pion masses ranging from 380MeV to 680MeV. We show that in this situation its performance is comparable to the recently proposed HMC variant with domain decomposition as preconditioner. We gi...
متن کاملNon-Hermitian Polynomial Hybrid Monte Carlo
In this thesis algorithmic improvements and variants for two-flavor lattice QCD simulations with dynamical fermions are studied using the O(a) improved Dirac-Wilson operator in the Schrödinger functional setup and employing a hybrid Monte Carlo-type (HMC) update. Both, the Hermitian and the Non-Hermitian operator are considered. For the Hermitian Dirac-Wilson operator we investigate the advanta...
متن کاملSolution of the Dirac equation in lattice QCD using a domain decomposition method
Efficient algorithms for the solution of partial differential equations on parallel computers are often based on domain decomposition methods. Schwarz preconditioners combined with standard Krylov space solvers are widely used in this context, and such a combination is shown here to perform very well in the case of the Wilson–Dirac equation in lattice QCD. In particular, with respect to even-od...
متن کامل